Study of Bs ^l^i. Decay in a Single Universal Extra Dimension 



o 

(N 



Ph. 



o 



Ying Li* Juan Hua, 
Department of Physics, Yantai University, Yantai 264-005, China 

January 21, 2013 



Abstract 

Utilizing form factors calculated within the light-cone sum rules, we have evaluated the decay branching ratios 
of Bs <j)Y and Sj — !> (j)£^£^ in a single universal extra dimension model (UED), which is viewed as one of the 
alternative theories beyond the standard model (SM). For the decay 6j — > (j)£^£^ , the dilepton invariant mass spectra, 
the forward-backward asymmetry, and double lepton polarization are also calculated. For each case, we compared the 
D ■ obtained results with predictions of the SM. In lower values of the compactification factor 1 /R, the only parameter in 

this model, we see the considerable discrepancy between the UED and SM models. However, when l/R increases, 
the results of UED tend to diminish and at 1/^ = 1000 GeV, two models have approximately the same predictions. 
Compared with data from CDF of 6, (j>fi^ fi^ , the l/R tends to be larger than 350 GeV. We also note that the 
fT^ ' zero crossing point of the forward-backward asymmetry is become smaller, which will be an important plat to prob 

' the contribution from the extra dimension model. The results obtained in this work will be very useful in searching 



rn 

' new physics beyond SM. Moreover, the order of magnitude for branching ratios shows a possibility to study these 

channels at the Large Hadron Collider (LHC), CDF and the future super-B factory. 



^ , 1 Introduction 

rS 
H 

. . .' As one kind of the flavor changing neutral current processes, the electroweak penguin decays b st^i appearing 
only at the loop level in the standard model (SM), are therefore sensitive to the fine structure of SM and to the possible 
new physics as well, and are expected to shed light on the existence of new physics before the possible new particles 
are produced at colliders. With the data from Belle and BaBar experiments, the exclusive processes B K^*''£^i^ 
have received great interest so that their theoretical calculation has been the subject of many investigations in the SM 
[1,2] and beyond [3, 4]. Along this line, the exclusive decays Bs (p£^£^ become also attractive since these decays 
are also induced by ^ — > s£'^£^, and could be measured at the running Tavatron, LHC and future super-B factories. 
Recently, the CDF collaboration had observed the rare semi-leptonic decay Bs (^pL^pL'^ [5] and the branching ratio 
is 

Br(B, ^ = [1.44±0.33(stat.)±0.46(syst.)] x 10"^. (1) 
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This exclusive process is quite worthy of intensive research and have attached much attention [6, 7]. When studying 
the semi-leptonic decays, dilepton invariant mass spectrum, the forward-backward asymmetry, and double lepton 
polarization are important observables to test SM and prob new physics, while the first two are mostly analyzed. 
Due to small mass of electron and muon, the invariant mass spectra and branching ratios are almost the same for 
electron and muon modes. Meanwhile, it is very difficult to measure the electron polarization, so we only consider 
Bs (/)jU+^",T+T" in this work. 

To make the theoretical predictions clearly, additional knowledge of decay form factors is needed, which is re- 
lated with the calculation of hadronic transition matrix elements, and can only be reliably calculated by using a 
non-perturbative QCD method. Fortunately, these form factors related to B, (j)£^£^ have been explored within the 
different methods, such as light cone sum rules (LCSRs) [8, 9, 10], perturbative QCD approach [11] and in different 
quark models: relativistic constitute quark model [12], constituent quark model [13], light front quark model [14]. 
Among them, the light-cone sum rules, which deal with form factors at small momentum region, is complementary 
to the lattice approach and has consistence with perturbative QCD and the heavy quark limit. We will adopt the form 
factors calculated by the LCSRs. 

Among the ideas proposed to extend the SM, a lot of attention has recently been devoted to models including extra 
dimensions [15, 16]. An interesting model is that proposed by Appelquist, Cheng and Dobrescu with the so-called 
universal extra dimensions (UED) [17], which means that all the SM fields may propagate in one or more compact 
extra dimensions. The compactification of the extra dimensions involves the appearance of an infinite discrete set of 
four dimensional fields which create the so-called KK particles. The simplest UED scenario is characterized by a 
single extra dimension. Compared to the SM, this model has one extra parameter called compactification radius, R. 
Hence, this model is a minimal extension of the SM in 4 + 1 dimensions with the extra dimension compactified to the 
orbifold /Z2 and the fifth coordinate, y running between and InR, and y — and y = nR are fixed points of the 
orbifold. The zero modes of fields propagating in the extra dimension correspond the SM particles. The masses of 
KK particles are related to compactification radius according to the relation =mQ + n^/R^ , with « = 1,2,.... One 
of the important property of the model is the conservation of KK parity that guarantees the absence of tree level KK 
contributions to low energy processes occurring at scales much smaller than the compactification scale. 

After the UED model being proposed, many attempts have been done to constraint the only parameter compactifi- 
cation radius R, for example, from Tevatron experiments the bound on the inverse compactification radius is found to 
be about l/R> 300GeV. The anomalous magnetic moment of muon and Z^bb vertex also lead to the same conclu- 
sion. Rare B transitions can also be used to constrain this scenario. Buras and collaborators [18] have investigated the 
impact of universal extra dimensions on the B[J s^B^d s mixing mass differences, on the CKM unitarity triangle and 
on inclusive b s decays for which they have computed the effective Hamiltonian. In particular, it was found that 
BR{B Xsj) allowed to constrain l/R > 250GeV, a bound updated by a more recent analysis to l/R > 600GeV at 
95% CL, or tol/R> 330GeV at 99% CL [19]. In this work, we will consider the 1 /R from 200 GeV up to 1000 GeV. 
In the past years, the UED model has been applied widely to calculate many observables related to the radiative and 
semileptonic decays of hadrons (see for example [4, 7, 20, 21, 22, 23]). 
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The aim of the paper is to find the effects of the KK modes on various observables related to the — > (pi'^i^ tran- 
sition, and these observables involve the dilepton invariant mass spectra, the forward-backward asymmetry, and double 
lepton polarization. We will compare the obtained results with the predictions of the standard model. In Ref. [7], R. 
Mohanta and A. K. Giri had calculated the processes (j)£^£^ and the Bs yi'^t' under the UED model, where 
they calculated the branching ratios and forward-backward asymmetries adding the long distance contribution. In this 
work, we will drop the contribution from the resonances, such as //i//, i//', and recalculate all observables. More- 
over, we will also calculate the lepton polarizations, which are always viewed as good places to prob the new physics 
contribution. In other words, this work can be regarded as supplementary of Ref.[7]. 

The outline of the paper is as follows. In section 2, after introducing the effective Hamiltonian responsible for 
the h st^ transition and form factors of Bj (j), we will present the formula of observable. In section 3, 
we numerically analyze the considered observables of B^ — > ^ji^fi^, Bs 0t+t^. This section also includes a 
comparison of the results obtained in UED model with that predicted by the SM. We will summarize this work at last. 



2 Effective Hamiltonian, Form Factors and Formula of Observable 

At quark level, the Bs (p£'^£^ transition proceed via FCNC transition of the b sl^i^ . Neglecting the doubly 
Cabibbo-suppressed contributions, the effective Hamiltonian governing b si'^i'^ transition is given by [24, 25] 

Jf^eff = y^V.hV* ^ C,(Ai)0,(Ai) , (2) 

v2 ,= 1 

where explicit expressions of (9/ could be found in Ref. [24], and the Wilson coefficients C,- can be calculated pertur- 
batively [26, 27, 28, 29]. Using the effective Hamiltonian, the free quark decay amplitude can be written as: 



GFa,mVtbV*s 



(3) 



Here q = p+ + p-, where p± are the four momenta of the leptons, respectively. Noted that ./#, although a free quark 
decay amplitude, contains part of long-distance effects from four-quark interactions, which usually are absorbed into a 
redefinition of the short distance Wilson coefficients. To be specific, we define the effective coefficient of the operator 
0<) as 

<f/-^ =C9 + Y{q^), (4) 

where Y{q^) stands for the above mentioned contribution from four-quark interaction. The corresponding operators 
and Y{q^) can refer to Ref. [30]. In this work, as mentioned before, we have not consider the contribution mainly due 
to J/\jf and xif' resonances in the decay chain Bs <j)£'^£^, which could be vetoed experimentally. 

The main source of the deviation of the UED model and SM predictions on the considered observables are from 
Wilson coefficients Cg^^ and Cio, which can be expressed in terms of the periodic functions, F{xi, 1 /R) with 
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X, = mj/M^ and being the top quark mass. Similar to the mass of the KK particles described in terms of the zero 
modes (« ~ 0) correspond to the ordinary particles of the SM and additional parts coming from the UED model, the 
functions, F{xt,\/R) are also written in terms of the corresponding SM functions, Fo(x;) and additional parts which 
are functions of the compactification factor, I /R, i.e.. 



F{x,,l/R) = Fq{x,) + ^ F„{x,,Xn), 



(5) 



111 11 

where Xn = — j- and »z„ = — . The Glashow-Illiopoulos-Maiani (GIM) mechanism guarantees the finiteness of the 
functions, F{xt, 1 /R) and satisfies the condition, F{x,, l/R) ^ Fo{x,), when R^Q. As far as l/Ris taken in the order 
of a few hundreds of GeV, the Wilson coefficients differ considerably from the SM values. For explicit expressions of 
the Wilson coefficients in UED model see [4, 18]. In Fig. 1, we plot the Wilson coefficients C, (/ = 7,8,9, 10) versus 
1 /R and find that the impact of the UED on the Cg is small. The suppression of jC? | and jCg |, that for l/R = 300 GeV 
amount to 82% and 66% relative to the SM values, respectively. For |Cio|, it can be enhanced by 16% for l/R = 
300 GeV, which does not renormalize under QCD. The UED contribution to four quark QCD penguin operators are 
also neglected in this work because of rather smaller Wilson coefficients. 
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Figure 1: The wilson coefficients = 4.8GeV) plotted versus 1 /R. The constant lines are the SM values. 



For the process — > (l)£^£ , the nonvanishing matrix elements are: 



{(P{k)\sY,,{l-75)bmp)) 



-ie*^i{mB, + m^)A\{q^) + i{2p - q) ^{e* ■ q) 



+iq^{e*-q)^^ Ai{q-')^Aoiq-' 



Mb, + niA 



(6) 



(0(^)|5-(7^v/(l + 75)^|£.(p)) 



1^l 



2 2 

mt -mi 



{2p~q)^ 



(7) 



where is the polarization vector of the meson. By means of equation of motion, one can obtain several relations 
between form factors 



A,{q^ 



Mb, + Md, 



Miq') 



2^ "«'~"^A2(^2) 



(8) 



2m^ 2m^ 

and Ao(0) ~ Ai,{Q),Ti{Q) = 72(0). All sighs are defined in such a way to render the form factors are real and positive. 
We will use the results calculated by using the technique of the light-cone QCD sum rule approach [8], the updated 
results and q^ dependence of the form factors can be found in Ref. [9]. In Ref.[7], the authors found that the uncertain- 
ties from the form factors can change slightly from their corresponding central values in the low s region, and highly 
suppressed in large s region, we will not consider the effect of these uncertainties here. The physical range in s = q^ 



extends from Smin — '^rnr to s 



I Lw J max 



Keeping the lepton mass and adopting the same convention and notation as [2], we find that the dilepton invariant 
mass spectrum forfi^ — > decay is given as 

dr Glal„m^ 



\V*Vrh\uiS)D . 



(9) 



where the function D is defined as: 

'-2 J{s)^ 
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1 

u(sf 



|B|2(A - + 8^2 (,f +2m|)) + \F\\X - ^ + 8^ (.'-4m|)) 



3 ' ' ' V 3 



'4mj{2 + 2ml-s) 



1 

2^ 



Re {BC* ) (A - ^ ) ( 1 - »4 - S) 



With S = s/m\, mi> ^ m/ /ms and = /niB, the kinematic variables are defined as 

A = l+m\,+f -2s-2m^^{\+s) 



-^2 -^2 

2^A \RtiFH*) - Rt(GH*){\ - ml)] + ^sX\H\^ . (10) 



u{s) = WA(l-4- 



(11) 
(12) 



Combined the effective coefficients and form factors, the auxiliary functions A ,B ,C ,E ,F ,G and H are refereed to 
Ref. [2]. According to the definition of the forward-backward asymmetry (FBA), it is straightforward to obtain the 
expression of the normalized FBA as: 

dAt 



^FB 



ds 



D = u{S)s[Re{BE*) +Re{AF*)] 



We define the three orthogonal unit vectors in the center mass frame of dilepton as 

eL = P+, = — — , eT=eNXeL, 

\PK^P+\ 



(13) 



(14) 



which are related to the spin of lepton by a Lorentz boost. Then, the decay width of the (j)£^£ decay for any 
spin direction ii of the lepton, where « is a unit vector in the dilepton center mass frame, can be written as: 

= ^(5)o[l + {PieL+PNeN+PTeT) ■ «] (15) 

where the subscript "0" denotes the unpolarized decay width. Pi and Pj are the longitudinal and transverse polar- 
ization asymmetries in the decay plane respectively, and P^ is the normal polarization asymmetry in the direction 
perpendicular to the decay plane. 

The lepton polarization asymmetry P, can be obtained by calculating 

_ dT{n^ei)/ds-dT{n = -ei)/ds 

^ dT{h = et)/dS+dT{n = -et) /dS ' ^ ' 

By a straightforward calculation, we get 



PlD = Jl-4^\—Re{AE^) + - ^Rc{BF^) 

V ^" i 3 

\J MBG^ +CF') + -^Re(CG+) , (17) 

+lm{FH^){l-ml~s)~lm{GH^)X] +2m^;«,[Im(B£^) +Im(AF' )] | , (18) 

P^D = 1 4sRe(AB^) + ^^'""^''^ \-Re{BF^) + (1 - ml)Re{BG^)+sRs{BH')] 

4VS \ mi 



-^[Re(CF^)-(l-m^)Re(CG^)-iRe(C//^)] J. (19) 



m 



3 Numerical Analysis 

In this section, we will examine the above mentioned physics observables and study their sensitivity to the compact- 
ification factor, l/R, which is the most important parameter in the single universal extra dimension model. Now, the 
parameters of the SM in our calculation are listed as follows: 

MB, = 5.36 GeV, = 4.8 GeV, m,. = 1.4 GeV, tn^ = 1.02 GeV, = 0.1057 GeV, 

m^^ 1.7769 GeV;7i, = 172.4 GeV, mw = 80.4 GeV, mz = 91.18 GeV, sin^ = 0.23, 

cce,n = av(mz) = 0.118, \V*V,b\ = 38.5 x 10-^ Tb, = 1.46 x lO-'^s. (20) 



It is interesting to calculate the branching ratio of B^ ^7 in the UED model firstly, since this channel will be 
measured easily among the B^ rare decays. Within the effective Hamiltonian, see equation (2), the decay rate of this 
channel is read: 

-z2 I m2 



r(B. ^ 07) = -^^yMmlm^ M - ;/ |cf ?|ri(0)p (21) 
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Figure 2: The branching ratio of 07 changes with l/R, and horizonal band corresponds to the SM prediction. 
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Figure 3: The variation of the dilepton invariant mass spectra of ^M^/i (left panel) and Bj 0t+t (right 
panel) with (in units of GeV^). The solid line corresponds to the SM, dotted line, dashed line and dot-dashed line 
are for 1 /R = 200 GeV,500 GeV, 1000 GeV respectively. 

In the SM, we found that the BR(Bs — s- 0/) = (2.4 ± 0.4) x 10^^, the major uncertainty is from the error of the form 
factor Ti (0). Considering the contribution of the UED, we present the branching ratio dependence on compactification 
parameter \ /R 'm Fig. 2. There are considerable discrepancies between the predictions of the UED and SM models for 
low values of the 1 /R. For l/R = 300 GeV, the ratio can reach to 1 .5^q3 x 10^^. Such a discrepancy at low values 
of 1/R can be a signal for the existence of extra dimensions. Accordingly, if the experiment can measure this channel 
well, we can constraint the range of 1 /R. 

In term of Eq. (10), we illustrate the dilepton invariant mass spectra in Fig. 3. By integrating the differential ratios 
over q^, in the SM, we obtain that: 

Br(B, ^ 0jU+Ai") = 1-8 X 10"^ 

Br(B,^0T+T") =2.4x10"^. (22) 

Considering the uncertainties in both theoretical side and experimental sides, the prediction of B^ (j)iJ.^ ji^ is con- 
sistent with CDF measurement Br(Z?s 0/i+/i") = (1.44±0.56) x 10^^ well, which is shown in Fig. 4. Note that 
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Figure 4: The branching ratio of Bs panel) and (right panel) changes with l/R, and 

horizonal solid lines corresponds to the SM prediction. In the left panel, the band area depicts the experimental data 
from CDF. 




Figure 6: The zero of the FBA changes with 1 /R, and horizonal line corresponds to the SM prediction. 
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the small branching ratio of tau mode is due to highly suppressed phase space. We find there is difference between our 
results and these from the Ref. [7], that is because we ignore the resonance contribution and slight different parameter 
space. Adding the the UED contribution, from the figures, one can see that the extra dimension effects can increase 
both the differential width and branching ratios for low values of \/R. As \/R increases, this difference tends to 
diminish so that for higher values of \/R{\/R> lOOOGeV), the predictions of UED become very close to the results 
of SM . Such discrepancy at low values of \/R can be considered as a signal for the existence of extra dimensions. 
Moreover, once l/R> SOOGeV, it is impossible to disentangle the extra dimension contribution clearly in both modes. 
As it is expected, the order of magnitudes of the branching ratios show a possibility to study such channels at the LHC. 

We show the lepton forward-backward asymmetries (FBA) for both channels in Fig. 5. It is shown that there 
is also considerable discrepancy between the predictions of the UED and SM models for low values of 1 /R. As 
l/R increases, this difference starts to diminish. It should be stressed that the hadronic uncertainties almost have no 
influence on this symmetry, so it can provide the best tool to prob the new physics contribution. For the muon mode, 
due to the destructive interference between the photo penguin and Z penguin, the FBA should have a zero point sq. 
Any deviation of the zero position sq from that of the SM would give us the clue of the new physics. In Fig. 5, one can 
find that sq decreases from the SM value. To make the situation clearer, the dependence of sq on the compactification 
parameter l/Ris depicted in Fig. 6, which shows that the zero position is pushed to a smaller value by decreasing 1 /R. 
Such a sensitivity indicates is particularly suited to constrain the R. Hence, with the enhancement of experimental 
precision and statistics in LHC and CDF, the measurements of FBA would provide more data and effectively pine the 
NP effects. 

Now, we turn to discuss the lepton polarization. In Fig. 7 and 8, we present the longitudinal and transverse 
polarization for Bs (pfi^fi^ and 0t+t^, respectively. For the normal polarization part, due to real Cio, it is 
the order of 10^^ and cannot be observed even in the designed super-B factory, then we will ignore this part in this 
work. From the Fig. 7 of Bs (pj^^jj,^, for small l/R, the deviation between UED model and SM is apparent. At the 
small value of momentum transfer, the \Pl\ is larger than the SM prediction; while at > 3.2GeV^, it will become 
smaller than that of SM. Thus, measurement of will discriminate between the different models. For the Pt part, the 
predictions from different models are almost the same. As B^ ^t+t^ as concerned, with smaller l/R = 200GeV, 
we can observe the deviation appears when > 16GeV^. Again in the Pt part the predictions are the same as SM for 
afl 1 /R. 

4 Summary 

Using form factors calculated within the light-cone sum rules, we presented the decay branching ratios of B^ ^ ^7 
and Bs — > ^£^£^ in a single universal extra dimension model (UED), which is view as one of the alternative theories 
beyond the standard model (SM) with only one new parameter For the decay Bs — > (pt^l^ , the dilepton invariant mass 
spectra, the forward-backward asymmetries, and double lepton polarizations are also calculated. For each case, we 
compared the obtained results with predictions of the SM. In lower values of the compactification factor 1/7?, we can 
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Figure 7: The dependence of longitudinal (left panel) and transverse(right panel) polarization of By — > ^ji^ji^ with 
(in units of GeV^). The solid line corresponds to the SM, dotted line, dashed line and dot-dashed line are for 
1 /R = 200 GeV,500 GeV, 1000 GeV respectively. 
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Figure 8: The dependence of longitudinal (left panel) and transverse(right panel) polarization of Bs T/i+/i^ with 
q-^ (in units of GeV^). The solid line corresponds to the SM, dotted line, dashed line and dot-dashed line are for 
1 /R = 200 GeV, 500 GeV, 1000 GeV respectively. 
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see the considerable discrepancy between the UED and SM models. However, when 1 /R increases, the results of UED 
tend to diminish and at l/R= 1000 GeV, two models have approximately the same predictions. For (pfi'^ji^, 
compared with data from CDF, the 1 /R tends to be larger than 350 GeV. We also note that the zero point of the 
forward-backward asymmetry become smaller, which will be an important plat to prob the contribution from the extra 
dimension. Moreover, the order of magnitude for branching ratios shows a possibility to study these channels at the 
Large Hadron Collider (LHC), CDF and the future super-B factory. 
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